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Abstract 



An aliTLOst-Riemannian structure on a surface is a generalized Riemannian structure whose 
local orthonormal frames are given by Lie bracket generating pairs of vector fields that can 
become coUinear. The distribution generated locally by orthonormal frames has maximal rank 
at almost every point of the surface, but in general it has rank 1 on a nonempty set which 
r^ ' is generically a smooth curve. In this paper we provide a short introduction to 2-dimensional 

"ti . almost-Riemannian geometry highlighting its novelties with respect to Riemannian geometry. 

Wc present some results that investigate topological, metric and geometric aspects of almost- 
Riemannian surfaces from a local and global point of view. 

'> ■ 1 Introduction 

0\ 

Q\^ ' The purpose of this paper is to present a generahzation of Riemannian geometry that naturally 

^D ■ arises in the framework of control theory. A Riemannian distance on a smooth surface M can be 

Cf^ I seen as the minimum-time function of an optimal control problem where admissible velocities are 

^^ ' vectors of norm one. The control problem can be written locally as 

^ q = uX{q) + vY{q), u^ + v^ < 1 , (1) 

^ I where {X, Y} is a local orthonormal frame. Almost-Riemannian structures (ARSs for short) gen- 

H ' eralize Riemannian ones by allowing X and Y to be collinear at some points. In this case the 

- - - corresponding Riemannian metric has singularities, but under generic conditions the distance is 

well-defined. For instance, if the two generators satisfy the Hormander conditioro (see for instance 
[6l dni [25l [30] ) , system ([1]) is completely controllable and the minimum-time function still defines 
a continuous distance on the surface. 

Let us denote by A(g) the linear span of the two vector fields of a local orthonormal frame at a 
point q. Where A(g) is 2-dimensional, the corresponding metric is Riemannian. Where A{q) is 1- 
dimensional, the corresponding Riemannian metric is not well defined, but thanks to the Hormander 
condition one can still define the Carnot-Caratheodory distance between two points, which happens 
to be finite and continuous. Generically, the singular set Z = {q € M \ dim(A(g)) = 1} is 
a 1-dimensional embedded submanifold and there are three types of points: Riemannian points, 



^The Hormander condition for a family T C Vec(M) of vector fields states that for every q £ M Lie^J^ — TqM. 



Grushin points where A{q) is l-dimensional and dini(A(g) + [A,A](g)) = 2 and tangency points 
where dini(A(g) + [A, A](g)) = 1 and the missing direction is obtained with one more Lie bracket. 
Generically, the following properties are satisfied: at Grushin points A(g) is transversal to Z, at 
tangency points A(g) is tangent to Z and tangency points are isolated. 

Almost-Riemannian structures on surfaces were introduced in the context of hypoelliptic oper- 
ators |221 I23j . They appeared in problems of population transfer in quantum systems [17^ [TSl [16] 
and they have applications to orbital transfer in space mechanics |12[lllj. 

The presence of a singular set enriches almost-Riemannian structures with several novelties 
with respect to the Riemannian case. The aim of this paper is to present and discuss some of these 
aspects, mainly from a geometric point of view. 

For instance, as it happens in sub-Riemannian geoemetry, spheres centered at points of the 
singular set are never smooth and the asymptotic of the Carnot-Caratheodory distance is highly 
non-isotropic, see Section 14.21 Moreover, cut loci accumulate at singular points, even in a non- 
smooth way, see [131 [18]. From the local point of view, the relations between curvature and conju- 
gate points change, as the presence of a singular set permits the conjugate locus to be nonempty 
even if the Gaussian curvature is negative, where it is defined. From the global point of view, the 
relations between curvature and topology of the surface change as well. Indeed, Gauss~Bonnet-type 
formulas for ARSs were obtained [S] [HI [8], where the role of the topology of the surface is instead 
played by the topology of the components in which the singular set split the surface and by the 
contributions at tangency points, see Section 15. li As a consequence, there exist ARSs on surfaces 
not necessarily parallelizable for which the integral of the curvature (suitably defined) vanishes, 
when in the Riemannian case the total curvature can vanish only on the torus. Also, the location 
of the singular set and of tangency points plays a fundamental role when studying the Lipschitz 
equivalence class of the Carnot-Caratheodory, see [19]. In particular, if in the Riemannian case 
all distances on the same surface are Lipschitz equivalent, in the almost-Riemannian case this is 
no longer true and the Lipschitz classification is finer than the differential one. Other interest- 
ing phenomena, not considered here, concerning the heat and the Schrodinger equation with the 
Laplace-Beltrami operator on an almost-Riemannian surfaces were studied in [20]. In that paper 
it was proven that the singular set acts as a barrier for the heat fiow and for a quantum particle, 
even if geodesies can pass through the singular set without singularities. 

The structure of the paper is the following. In Section [2] we recall the precise definition of almost- 
Riemannian structure on a surface and fix some notations. In Section [3] we provide examples: the 
Grushin plane, for which we compute the optimal synthesis, and an example on the 2-sphere. Sec- 
tion [J] deals with local results. First in Section [4.11 we study local orthonormal frames and state a 
result that classifies the types of points. Then, in Section [12] we focus on local properties around 
tangency points such as the asymptotic of the distance as well as the cut and conjugate locus. 
Section [4.31 discusses functional invariants, i.e., functions that allow to recognize locally isometric 
structures. In Section [5] we address global problems. First we present in Section [5.11 several gener- 
alizations of the Gauss-Bonnet formula in the almost-Riemannian context. (For generalizations of 
Gauss-Bonnet formula in related contexts, see also [5| I27[[28] .) In Section [5.21 a characterization of 
ARSs using the topology of the vector bundle defining the structure is stated. Finally, in Section [53] 
we analyse almost-Riemannian surfaces from a metric point of view presenting a classification result 
with respect to Lipschitz equivalence. We conclude in Section [6] by some open questions. 



2 Basic definitions and notations 

Let M be a smooth connected surface without boundary and denote by Vec(M) the module of 
smooth vector fields on M. Throughout the paper, unless specified, objects are smooth, i.e., of 
class C~. 

Definition 1 An almost-Riemannian structure (ARS for short) on a surface M is a triple {E, f, (•, •)), 
where (i) E is a Euclidean bundle of rank two over M (i. e. a vector bundle whose fibre is equipped 
with a smoothly-varying scalar product {■, ■)q); (ii) f : E ^ TM is a morphism of vector bundles 
such that f{Eq) C TqM ; (Hi) for every q £ M LicqA = TqM, where 

A = {f o cr I cr section of E} . 

ARSs on surfaces can be seen as a first generalization of Riemannian structures towards sub- 
Riemannian ones. Indeed, recall that a (constant rank) sub- Riemannian structure on a manifold 
N is given by (D,g), where D C TN is a sub-bundle of rank k < dimA^ such that Lie^-D = TqN 
for every q £ N and g' is a Riemannian metric on D. When diniA^ = 2, the only possibility for a 
sub-bundle to be Lie bracket generating is that D = TN, i.e., there do not exist constant rank sub- 
Riemannian structures. Hence, in this case one needs to consider rank- varying modules A C VecA^ 
instead of sub-bundles D C TN. 

Let S = {E,f, {■,■)) be an ARS on a surface M. The Lie bracket generating assumption in 
Definition [H jointly with the fact that diniM = 2, implies that for each vector field X G A there 
exists a unique section a of E such that X = f o a. 

Denote by A{q) the linear subspace {V{q) \ V £ A} = f{Eq) C TqM. The set 

Z = {q £ M \ dim A{q) < 2} 

is called singular set and coincides with the set of points q £ M at which fl^;,^ is not injective, i.e., 
dimA(g) = 1. 

A property (P) defined for 2- ARSs is said to be generic if for every rank-2 vector bundle E over 
M, (P) holds for every f in an open and dense subset of the set of morphisms of vector bundles 
from E to TM inducing the identity on M, endowed with the C°°- Whitney topology. 

The Euclidean structure (•,•) on E induces a symmetric positive definite bilinear form G{-,-) 
on the submodule A as G{V,W) = ((Ty,(Tiy), where av,crw are the unique sections such that 
V = fo av,W = f o aw- At points q where fl^;^ is an isomorphism G acts as a tensor, i.e., 
G{V,W)\q depends only on V{q),W{q). This is no longer true at points belonging to Z, which is 
generically a smooth embedded submanifold of dimension 1. By definition, an ARS is Riemannian 
if and only if f is an isomorphism of vector bundles or, equivalently, the singular set is empty. 

If ((Ti,(T2) is an orthonormal frame for (•, •) on an open subset U of M, an orthonormal frame 
for 5 on f7 is given by (f o ai, f o (T2). 

For every q £ M and every v £ A{q) define Gq{v) = ini{{u,u)q \ u £ Eq,f{u) = v}. If q ^ Z, 
then Gq{v) = G{V,V)\q, for any vector field V such that V{q) = v. li q £ Z then we have the 
inequality 

G{V,V)\q>Gq{V{q)). 

An absolutely continuous curve 7 : [0, T] — )• M is admissible for S if there exists a measurable 
essentially bounded function [0, T] 9 i 1— )■ u{t) £ E^(^i-j such that 7(t) = f{u{t)) for almost every 



^(7) = / ^G,it){i{t)) dt. 



t E [0,T]. Given an admissible curve 7 : [0,T] — )• M, the length of 'j is 

r-T 
'0 

The Carnot-Caratheodory distance (or almost-Riemannian distance) on M associated with S is 
defined as 

d{qo,qi) = inf{£(7) | 7(0) = qo,j{T) = qi,j admissible}. 

The finiteness and the continuity of d{-,-) with respect to the topology of M are guaranteed by 
the Lie bracket generating assumption (see [9, Theorem 5.2]). The Carnot-Caratheodory distance 
endows M with the structure of metric space compatible with the topology of M as differential 
manifold. 

Locally, the problem of finding a curve realizing the distance between two fixed points qo,qi G M 
is naturally formulated as the distributional optimal control problem 

q = uiFi{q) + U2F2{q) , Ui G M, / ^ul{t) + ul{t) dt -^ min, g(0) = qo, q{T) = qi, 

where Fi,F2 is a local orthonormal frame for S. 

A geodesic for S is an admissible curve 7 : [0, T] — )• M, such that GryU\{'y{t)) is constant and 
for every sufficiently small interval [ti,t2] C [0, T], 7|[ii,i2] ^^ ^ minimizer oi i. A geodesic for which 
G^(i)(7(t)) is (constantly) equal to one is said to be parameterized by arclength. 

Although the metric tensor is singular at points of Z, geodesies are well-defined and smooth. 
This can be proved by using classical methods of optimal control theory. The Pontryagin Maximum 
Principle [29] provides a direct method to find geodesies for an ARS, as the Hormander condition 
ensures the absence of abnormal extremals [3, Proposition 1]. Indeed, it implies that an admissible 
curve parameterized by arclength is a geodesic if and only if it is the projection on Af of a solution 
of the Hamiltonian system corresponding to the Hamiltonian 

H{q,p) = \{ip-F,iq))^ + ip-F2iq)f), q G U, p G T*U. (2) 

lying on the level set H = 1/2, where Fi,F2 is a local orthonormal frame for the structure on an 
open set U. Notice that H is well defined on the whole T*M, since formula ([2]) does not depend 
on the choice of the orthonormal frame. When looking for a geodesic 7 realizing the distance 
from a submanifold N (possibly of dimension zero), one should add the transversality condition 
p{0) ± T^(o)iV. 

The cut locus Cut^ from A'^ is the set of points p for which there exists a geodesic realizing 
the distance between A^ and p that loses optimality after p. It is well known (see for instance |j] 
for a proof in the three-dimensional contact case) that, when there are no abnormal extremals, if 
p £ Cut/vT then one of the following two possibilities happen: (i) p is reached optimally by more 
than one geodesic; (ii) p belongs to the first conjugate locus from N defined as follows. To simplify 
the notation, assume that all geodesies are defined on [0,oo[. Define 

C° = {A = {q,p) e T*M \qeN, H{q,p) = 1/2, p ± T,N] 



and 

exp: C° X [0,oo[ -^ M 



tH ^ 



{X,t) ^ 7r(e'^A) 

where vr is the canonical projection {q,p) — )> q and H is the Hamiltonian vector field corresponding 
to H. The first conjugate time for the geodesic exp(A, •) is 

tconjW = min{t > 0, (A,t) is a critical point of exp}. 

and the first conjugate locus from N is {exp(A, tconj(A)) | A G C^}. 

3 Examples 

3.1 The Grushin plane 

To have a better understanding of the topic, let us study the simplest case of genuinely almost- 
Riemannian surface. 

Consider the ARS on M^ ^^ere E = R'^ x M^ f((x,y), (a, 6)) = {{x,y),{a,xb)), and (•,•) is 
the canonical Euclidean structure on M^. In this case a global orthonormal frame is given by 
Fi{x,y) = dx, F2{x,y) = xdy and the singular set is indeed nonempty, being equal to the y-axis. 
This ARS is called Grushin plane, named after V.V. Grushin who studied in [23j analytic properties 
of the operator d^ + x'^dy and of its multidimensional generalizations (see also |22j). 

The bilinear form G{-,-) in coordinates (x,y) reads 

1 

^ 

and the Gaussian cuvature is K{x, y) = —\. Note that for every point outside the singular set the 
curvature is negative and it diverges to — oo as the point approaches Z. 

The Grushin plane provides a very illustrative example as geodesies can be computed explicitly. 
Indeed, by the Pontryagin Maximum Principle, geodesies are projection on M^ of solutions of the 
Hamiltonian system associated with 



H{x,y,px,Py) = -{pI+x^pI). 



The system is 



• 

x 


= 


Px 


y 


= 


X^Py 


Px 


= 


-xpI 


. Py 


= 






The last equation implies that along a geodesic Py{t) = Py{0) = a G M. For a = 0, we obtain 

x{t) = x{0) + px{0)t, y{t) = 0, Px{t)=Px{0), 



that is, horizontal half-hnes are geodesies. For a ^ 0, the first and the third equation are the 
equations for an harmonic oscillator whose solution is 

x(0) 

x{t) = x(0)cos(at) + — ^sin(at) (3) 

a 

Px{t) = — a3;(0) sin(at) + x(0) cos(at). 
Using the normalization H = 1/2, the initial covector {px{0),a) satisfies 

H{x{0),y{0),p^{0),pym = l{pl{0)+xH0W) = \, 

whence 

i(0) = Px(0) = ±Vl-a;(0)2a2. 

Integrating the equation for y(-) we get 

y{t) = y^,^^mm^aMof+m\^ 

a2x(0)2 - i(0)2 ^ ^ x(0)i(0) , , 

+ ^-^ ^^ sin(2at L^_^ cos 2at . (4) 

2a 2a 

Taking (0,0) as starting point, the normalization condition H = 1/2 implies Px{^) = 1 or 
Pa:(0) = —1- When ^^^(O) = 1, respectively Px{^) = —1, the geodesic enters the region {{x,y) \ x > 
0}, respectively {{x,y) \ x <0}. Hence, geodesies starting at (0,0) have the form (x"''(t, a),y(t,a)) 
or (x~(t, a),y(t,a)), where 

1 11 

X (t,a) = lb— sin(at), y{t,a) = — t ;r sin(2at) a / 0, (5) 

a 2a 4a^ 

x^{t,0) = ±t, y(t,0) = 0, a = 0. (6) 

Geodesies with a > 0, respectively a < 0, are contained in the half-plane {{x,y) \ y > 0}, respec- 
tively {{x,y) I y < 0}. In Figured^ geodesies for some values of a are portrayed, while Figure [TJd 
illustrates the set of points reached in time t = 1. Notice that this set is non-smooth. Also, the 
spherqj centered at (0, 0) of radius r = 1 is non-smooth. In contrast with what would happen 
in Riemannian geometry, this is the case for every positive radius, as it happens in constant-rank 
sub-Riemannian geometry. However, this is a consequence of the fact that the initial condition 
belongs to Z. 

Let us compute the cut locus from (0,0). The geodesies with a = never lose optimality. 
Due to the symmetries of the problem, it is easy to see that the two geodesies {x^{t,a),y{t,a)), 
(x~(t, a),y{t, a)) are optimal until they intersect at time t = 7r/\a\. As a consequence the cut locus 
from the origin is the set Cut(o,o) = {(0, a) | a € M\{0}}. Note that Cut(o,o) accumulates at (0, 0). 
This is due to the fact that (0, 0) & Z and represents another difference with the Riemannian case, 
where Cutp is always separated from p. 



^The sphere centered at (xq, yo) of radius r is defined as the boundary of the set {{x, y) \ d{{x, y), [xq, yo)) < r}, 
where d is the Carnot-Caratheodory distance. In general, for ARSs this set contains {{x,y) \ d{(x,y), (a;o,j/o)) = r} 
as a proper subset. 





Figure 1: Analysis of the Grushin plane at (0,0). Figure a shows some geodesies starting at (0,0) 
corresponding to a = ±1, ±2, ±3. Figure b portraits the set of points whose distance from (0, 0) is 
equal to 1 



To compute the conjugate locus from the origin, we find critical points of the map exp defined 
in Section [2j Using ([5j) one can check that, for a ^ 0, the first conjugate time on a geodesic 
(x^(t, a),y(t, a)) coincides with the first positive root of the equation 

atcos(af) — sin(at) = 0. 

Therefore the conjugate locus from (0, 0) is the union of two parabolas 

(x, y) \\y\ = ^ ( ^— - -) I \ {(0, 0)}, 

2 \ cos T sm T T J ] 

where r is the first positive number such that tanr = r. Note that since (0,0) G Z, the conjugate 
locus accumulates at (0,0). In particular, even if the curvature is negative for every point outside 
Z the conjugate locus from (0, 0) is non-empty. However, every geodesic reaches its conjugate time 
after crossing the singular set, see Figure [3ti. 

One may infer that the existence of conjugate points depends on (0, 0) belonging to the singular 
set. However, this is not the case. Indeed, let us consider (—1,0) as starting point. Let us 
parameterize the initial covector as Px{Q) = ±Vl — a^, Pj/(0) = a. When a = the solutions of the 
Hamiltonian system are x^(t,0) 
are 



-1 ± t, y{t,0) = 0. When a / 0, using ([SD, ^, the solutions 



X {t,a) 
y^{t,a) 



( cos (at) ± Vl — a^ sin(at) 



VT^^ 



{cos{atf - 1) + 



2at - sin(2at) + 2a'^ sin(2at) 



4a2 



^), (x+, -y+) (when ^^,(0) = 
^'^ As one expects, since the 



For a 7^ 0, geodesies starting at (—1,0) are curves of the form (x" 

Vl — a^) and of the form {x^ ,y^), (x~, —y~) (when PxiO) = — VT— a 

metric is Riemannian at (—1,0), for r < 1 the sphere centered at (—1,0) of radius r is smooth. 

This is no longer true for r > d((— 1, 0), (0, 0)) = 1, see Figure [21 

To compute the cut locus from (—1,0), note first that the horizontal geodesies {x^{t,0),0) 
never lose optimality. On the other hand, when a 7^ 0, one can check that the first time at which 
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-40^ 



Figure 2: The set of points {{x,y) \ d{{x,y), (—1,0)) = r}, for r = 1, 10. 



conjugate 
locus 




conjugate 
locus 




Figure 3: Geodesies and upper part of first conjugate locus from a point in Z (Figure a) and from 
a point outside Z (Figure b) 



a geodesic {x~^{t,a),y^{t,a)) intersects another geodesic, namely {x~{t,a),y^{t,a)), is t = vr/|a|. 
Similarly, the first intersection of (x+(t, a), — y"*"(t, a)) happens for t = vr/|a| with the geodesic 
(x~(t, a), —y~{t, a)). Hence, the cut locus from (—1, 0) is the union of two half-lines 



Cut 



(-1,0) 



{(l,a) I a G 



-oo, 



-V2)U(7r/2,+oo)}. 



As concerns the conjugate locus, computing critical points of exp it is easy to show that for 
each a 7^ every geodesic has a positive conjugate time. Hence, even if the curvature is negative 
and (—1,0) ^ Z, there exist conjugate points to (—1,0), see Figure [3)3. Notice that each geodesic 
reaches the conjugate point after crossing at least one time the singular set. 



3.2 An example on the 2-sphere 

Consider another example on a compact surface. Let M = S"^ C M'^. Then every pair of vector 
fields on S'^ are linearly dependent on a nonempty set. Hence, if one wants to study a metric 
structure defined globally by a pair of vector fields on the 2-sphere, then one needs to consider an 
almost-Riemannian structure. 

Metric structures defined globally by a pair of vector fields on S^ arise naturally in the context 
of quantum control (see \n\ 115]). Indeed, consider the ARS on S'^ where E is the trivial bundle 
of rank two over S'^ and the image under f of a global orthonormal frame for (•,•) on E is the 
pair X{x, y, z) = (y, —x, 0), Y{x, y, z) = (0, z, —y). Then the two generators are linearly dependent 



Integral Curves of X 




Figure 4: Alinost-Rieniannian structure on the 2-sphere 



on the mtersection of the sphere with the plane {{x,y,z) | y = 0} (see Figure H]). In this model, 
the sphere represents a suitable state space reduction of a three-level quantum system and the 
orthonormal generators X and Y are the infinitesimal rotations along two orthogonal axes, modeling 
the action on the system of two lasers in the rotating wave approximation. 

Note that, as for the Grushin plane, the distribution is transversal to the singular set at each 
point. Indeed, this ARS is the compact correspondent to the Grushin plane. 

4 Local results 

4.1 Local representations 

The first important work studying general properties of ARSs is [3j where the authors provide the 
characterization of generic ARSs by means of local representations. 

Definition 2 A local representation of an ARS S at a point q G M is a pair of vector fields {X, Y) 
on M? such that there exist: (i) a neighborhood U of q in M, a neighborhood V of (0, 0) in M? and 
a diffeomorphism ip : U — 
around q, such that (p^Fi 



V such that (p{q) = (0,0); (ii) a local orthonormal frame (i<i,i<2) of S 
= X, ip^:F2 = Y, where 99^, denotes the push-forward. 



Under generic assumptions, it turns out that one can always construct a local representation (X, Y) 
where the first vector field is rectified and the second one has a simple form. The main assumption 
to obtain such result is the following. Set Ai = A and Afc+i = A^ + [A, A^], i.e., A^ is the module 
spanned by Lie brackets of length less than k between elements in A. We say that S satisfies 
condition (HO) if the following properties hold: 

(HO) (i) Z is an embedded one-dimensional submanifold of M; 

(ii) the points q G M where A2((?) is one-dimensional are isolated; 
(iii) A3(g) = TqM for every q £ M. 

A simple transversality argument allows to show that property (HO) is generic for 2- ARSs. 



9 



Remark 1 Throughout the paper, unless specified, we always deal with ARSs satisfying (HO). 

Theorem 1 ([3j) Given an ARS S on M, for every point q €z M there exist a local representation 
{X,Y) of S at q such that {X,Y) has one of the forms 

(Fl) X{x, y) = (1,0), y(x, y) = (0, e^^^'^^), 

(F2) X{x,y) = (1,0), Y{x,y) = (0,xe'^(^'^)), 

(F3) X{x, y) = (1,0), Y{x, y) = (0, {y - x^^{x))e^^^'y^), 

where (p, ^ and ip are smooth real-valued functions such that 0(0, y) = and V'(O) 7^ 0. 

Thanks to Theorem [TJ for a point q € M there are three possibihties. First, if A(g) = TqM 
then q is said to be a Riemannian point and 5 is locally described by (Fl). Second, if A(g) is 
one-dimensional and A2{q) = TqM then we call q a Grushin point and the local description (F2) 
applies. At Grushin points ^{q) is transversal to Z and the Lie bracket between the two elements 
of a local orthonormal frame is sufficient to span the tangent plane TqM. Third, if A(g) = ^2{<l) 
has dimension one and A3(g) = TqM then we say that q\s a. tangency point and S can be described 
near q by (F3). At tangency points the subspace /S.{q) is tangent to Z and the missing direction is 
obtained with a Lie bracket of length two between the two elements of a local orthonormal frame. 
We also set 

T = {q ^ Z \ q tangency point of S}. 

Note that condition (ii) in assumption (HO) ensures that tangency points are isolated, i.e., T is a 
discrete set. 

The idea behind the proof of Theorem [1] is that, since A(g) is at least one-dimensional, given 
a local orthonormal frame for 5, one vector field can always be rectified. Then, to deduce the 
form of the other vector field one needs to construct a suitable coordinate system. To do this, 
the authors provide a procedure that allows to build a local representation starting from a smooth 
parameterized curve c(-) passing through the base point q and transversal to the distribution at 
each point. The transversality of c(-) to the distribution implies that the Carnot-Caratheodory 
distance from the support of c(-) is smooth on a neighborhood of q (see pi Lemma 1]). Given a 
point p near q, the first coordinate of p is, by definition, the distance between p and the chosen 
curve, with a suitable choice of sign. The second coordinate of p is the parameter corresponding to 
the point (on the chosen curve) that realizes the distance between p and the curve (see Figure [5]). 
Using the inverse of the diffeomorphism defining this coordinate system, one can build a vector 
field belonging to A of norm one and then complete it to a local orthonormal frame. The local 
orthonormal frame {X, Y) obtained through this procedure has the form X = (1, 0), Y = f{x, y)Y, 
where / is a smooth function. 

4.2 Local analysis at tangency points 

ARSs are characterized by the presence of a singular set, which includes Grushin and tangency 
points. These two types of points are essentially different from each other. 

At Grushin points, the local situation is described by the Grushin plane (see Section [3]) which 
is the nilpotent approximatiorl^ of any ARS at a Grushin point. The fact that the distribution 



■^The nilpotent approximation of a sub- Riemannian structure at a point is the analogue of the tangent space of a 
manifold. For the precise definition and analysis of the subject see for instance [6] |10ll30] . 
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sign =: —%sign. ^ 1 




sign — —1 sign ^ 1 



(a,a) [s.a) 



9- (0,0) 



(-s',a') (0,q:') 



c(-): curve transversal to the distributic 



Figure 5: The construction of coordinates starting from a parameterized curve c :] — e, £[—7- M. We 
denote by 7^ the geodesic starting at c{a), parameterized by arclength, entering the region where 
sign = =bl and such that d(7a (s),c(] — e,e[)) = s. As the distribution is transversal to c(-) at each 
point, the distance from c(] — e, e[) is smooth. 



is transversal to the singular set in a neighborhood of such points implies that the behaviour of 
the almost-Riemannian distance from Z is similar on the two sides of Z, as we point out in the 
discussion following Theorem [3] (see Section I5.1.ip . 

At tangency points, the situation is more complicated due to the fact that the asymptotic of 
the distance from the singular set is different from the two sides of the singular set. To see this, 
let us consider an example. Take the ARS on M? having X{x,y) = (1,0), Y{x,y) = (0, y — x^) 
as a global orthonormal frame. In this case, the singular set is the parabola {(x,y) \ y — x'^ = 0} 
(see Figure [6|) and (0, 0) is a tangency point. Define M^ as the set of points p £ M? such that 
d(p, Z) > e, where d is the almost-Riemannian distance. Then M^ is split in two connected 



components M> = M^ D {{x, y) \ y - x^ > 0} and M< = M^ H {{x, y) \ y 



< 0} contained on 



the two different sides of Z. It turns out that for small values of e both parts of the boundary of 
Me are non-smooth in a neighborhood of (0, 0) but have a corner. Also, the order at which they 
approach Z as e goes to zero is different. Indeed, taking another distance d on M^ associated with 
any Riemannian structurqj, one can check that d{M^,Z) = O(e^), whereas d{M^,Z) = 0{e^). 

The asymmetric behaviour of the almost-Riemannian distance distinguishes tangency points 
from Grushin points and has various consequences. It affects the shape of the cut locus from the 
tangency point, as well as the shape of the cut locus from Z near the tangency point. Moreover, 
together with the asymptotic of the Gaussian curvature, it is to be taken account of when proposing 
a notion of integrability with respect to an ARS (see Section [5. 1.2p . Let us analyse singular loci 
around tangency points. 

Consider an ARS in a neighborhood of a tangency point. Thanks to Theorem[Tl this is equivalent 
to taking the ARS on E? whose orthonormal frame is 

X{x,y) = {l,0), Y{x,y) = {0,{y-x^^{x))e^'^-'y^), (7) 

where ip, ^ are smooth functions depending on the structure and "0(0) 7^ 0- It is easy to see that 



'For example choose d as the Euclidean distance on " 
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Figure 6: The singular set (thick hue) and the boundary of M^ (bold line) in a neighborhood of a 
tangency point, for the structure on M^ whose orthonormal frame is (1,0), (0,y — x^) and e = 0.07 



the nilpotent approximation of d?]) at (0, 0) is the ARS defined by the orthonormal frame 

X(x, y) = (1,0), Y{x, y) = (0, -^x^), (8) 

where 7 = ■0(0)6^''^''^^ The singular set of this structure is the y-axis and at each singular point 
dimA(0,y) = 1. Therefore, this structure does not satisfy condition (ii) in (HO). The optimal 
synthesis for this ARS was computed explicitly in [21 [H] in terms of Jacobi elliptic functions. Even 
though such synthesis does not reflect qualitative properties of the one for the generic case ([7D , it 
can be used to compute the jet of the exponential map for the ARS ([7]). The development at (0, 0) 
of the orthonormal frame in ([7]) truncated at order zerqfl is 



X 



__,y = ^[y-x'^ 



e'x^^ 



d_ 

dy' 



(9) 



where e' = ip'{^) + V'(0)'^x(0, 0). The following proposition computes the exponential map at the 
origin for the ARS defined in ([9]). It happens that higher order terms in the expansion of the 
elements of the orthonormal frame in ([7]) do not affect the estimation of the exponential map and, 
consequently, the order zero is sufficient to describe the cut and conjugate loci from the tangency 
point, at least qualitatively. 

Proposition 1 ([13j) Consider the ARS on M? defined by the orthonormal frame given in ^. 
The solution of the Hamiltonian system, associated with 



H{x,y,px,Py) = ^{pl + l^{y-x^ 



e'x'fvl) 



with initial condition {x,y,Px,Py)\t=o = (0, 0,ibl,a) with 



a\ ~ 



+00 can be expanded as 



where rj 



x{t,ri) = ri x^{t/r]) +rf'x^{t/T]) + o{Tf), 
y{t,r,) = ry3yO(f/r?)+r?V(i/^)+o(^'), 

, {x^,y^,PxjPy) is the extremal of the nilpotent approximatior^ with initial condition 



^The coordinate functions {x,y) have weights (1,3), see 1101 . 

®that is, (a;", y'^,PxiPy) is a solution of the Hamiltonian system associated with ff" = l/2(p^ + x'^Py), see (|8]) 
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{x'^,y\pl,pl)\t=o = (0,0,±l,sign(a)), and 



x^ = pI, 



yi = 72(pi(x0)4 + 4p0(2;0)3xi-2p0((x0)2y0-e'(x0)5)), 

pI = -72(4p0pi(x0)3 + 6(p0)2(x0)22;i - {p^y)^2x'^y'^ - 5e'{xy)) 

pi, = TPy X , 

with initial condition {x^,y^,pl:,Py)\t=o = (0,0,0,0). 

Proposition [T] was used to estimate the conjugate locus from (0,0) for the ARS on M? defined 
by the orthonormal frame d?]), see [131 Proposition 5]. There exists a constant a ^ such that the 
conjugate locus from (0, 0) accumulates at (0, 0) as the set 

{{x, y)\y = ax^} U {{x, y) \ y = -ax^} \ {(0, 0)}. 

The shape of the cut locus from a tangency point (see Figure [7]) is described by the following 
result. 

Proposition 2 ([13j) Let S be an ARS on M and q € M be a tangency point such that there 
exists a local representation of the type (F3) for S at q with the property 

V^'(o) + 7/.(o)5,.e(o,o)/o. 

Then the cut locus from the tangency point accumulates at q as an asymmetric cusp whose branches 
are separated locally by Z. In the coordinate system where the chosen local representation is (F3), 
the cut locus accumulates as the set 

{{x, y) I y > 0, 2/2 - aix^ = 0} U {{x, y)\y<0,y'^ - a2X^ = 0}, 

with ai = Cic"^^^^'^' / {t/j' (0) + tlj{0)dx^{0,0))^ , the constants ci being nonzero and independent on the 
structure. 

Note that both the conjugate locus and Cut(o,o) accumulates at (0, 0) with tangent direction parallel 
to the distribution. This is no longer true for Cut^ in a neighborhood of a tangency point. A 
description of such locus is given by the following theorem, see also Figure [71 

Theorem 2 ( |18| ) Let S be an ARS on M and q £ M be a tangency point such that there exists 
a local representation of the type (F3) for S at q with the property 

e' = /(o) + v(o)a,e(o,o)/o. 

Then the cut locus from the singular set Z accumulates at q as the union of two curves locally 
separated by Z. One of them is contained in the set {y > x'^il){x)^, takes the form 

{(-l/2'0'(O)t + o(t), t + o{t)) I t > 0}, 

and accumulates at q transversally to the distribution. The other one is contained in the set {y < 
x'^'ip{x)} and takes the form 

{{e' ujt^ + o{t^), -t^ + o{t^)) I t > 0}, 

where u ^ is a constant depending on the structure. This part of the cut locus accumulates at q 
with tangent direction at q belonging to the distribution. 
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singular set 




cut locus- 




the cut locus from the tangency point 



the cut locus from the singular set 





singular set 



crest of the curvature 
■ transversal to the 
distribution 



crests and valleys of the curvature 



curves through the tangency point 



Figure 7: The singular set (dotted line), the cut locus from a tangency point (semidashed line), the 
cut locus from the singular set (dashed line), and the set of crests and valleys of K (solid lines) for 
the ARS with orthonormal frame Fi = -^, F2 = {y — x^ — x^)-^. In this case there are three curves 
in the set of crests and valleys of the curvature, only one of which is transversal to the distribution. 



4.3 Functional invariants 



In a recent paper [18] the authors address the problem of finding normal forms for ARSs that 
are completely reduced, in the sense that they depends only on the ARS and not on its local 
representation. This consists in finding a canonical choice for a local coordinate system and for a 
local orthonormal frame, i.e., two vector fields X and Y defined in a neighborhood of the origin on 
M?. Once a canonical choice of {X, Y) is provided, the two components of X and the two components 
of Y are functional invariants of the structure, in the sense that locally isometric structures have 
the same components. Moreover, they permits to recognize locally isometric structures: if two 
structures have the same invariants in a neighborhood of a point, then they are locally isometric. 
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Notice that the problem of finding a set of invariants that determines the structure up to local 
isometries is not completely trivial even in the simplest case of Riemannian points. See for instance 
the discussion in [Jl [26]. Indeed even if one is able to fix canonically a coordinate system, the 
Gaussian curvature in that coordinate system is an invariant, but there are non- locally isometric 
structures having the same curvature. 

A first step in finding normal forms is Theorem [TJ However, the local representations cor- 
responding to Riemannian and tangency points are not completely reduced. Indeed, there exist 
changes of coordinates and rotations of the frame for which an orthonormal basis has the same 
expression as in (Fl) (respectively (F3)), but with a different function </> (respectively with dif- 
ferent functions ip and ^). Recall that to build the local representations in Theorem [T] a specific 
procedure based on the choice of a smooth parameterized curve transversal to the distribution was 
used. If the parameterized curve used in this construction can be canonically built, then one gets 
a local representation of the form X = (1,0), Y = (0, /) which cannot be further reduced. Hence, 
/ is a functional invariant that completely determines the structure up to local isometries. In |18j . 
a canonical choice of the parameterized curve at each point is provided and the properties of the 
functional invariant / are studied. 

For Riemannian points, a canonical parameterized curve transversal to the distribution can be 
easily identified, at least at points where the gradient of the Gaussian curvature is non-zero: one 
can use the level set of the curvature passing from the point, parameterized by arclength. For 
points where the gradient of the curvature vanishes, under additional generic conditions, a smooth 
parameterized canonical curve passing through the point can be selected among the crests and 
valleys of the curvature, see [IHl Sections 4.1.3, 4.1.4]. 

For Grushin points, a canonical curve transversal to the distribution is the set Z. This curve 
has also a natural parameterization and was used in [3] to get the local representation (F2) that, 
as a consequence, cannot be further reduced. 

To obtain the local expression (F3), the choice of the smooth parameterized curve was arbitrary 
and not canonical. The analysis in [18] is aimed at finding a canonical one and, as a consequence, 
obtain a functional invariant at a tangency point that completely determines the structure. The 
most natural candidate for such a curve is the cut locus from the tangency point. Neverthless, by 
Proposition [2] this is not a good choice, as in general the cut locus from the point is not smooth 
but has an asymmetric cusp (see Figure [7]) . Another possible candidate is the cut locus from the 
singular set in a neighborhood of the tangency point, but Theorem [2] states that the cut locus from 
Z is non-smooth in a neighborhood of a tangency point, see Figure [Tj A third possibility is to 
look for curves which are crests or valleys of the Gaussian curvature and intersect transversally the 
singular set at a tangency point. The main result in [18] consists in the proof of the existence of 
such a curve. Moreover, this curve admits a canonical regular parameterization. 

An example of a crest of the curvature at a tangency point is shown in Figure [8] 

5 Global results 

In this section we address global questions for almost-Riemannian surfaces. The general context in 
which all the results can be stated is the one of totally oriented structures. 

Definition 3 An ARS is said to be oriented if E is oriented as vector bundle. We say that an ARS 
is totally oriented if both E and M are oriented. For a totally oriented ARS, M is split into two 
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for the almost-Riemannian structure 



Figure 8: The graph of the Gaussian curvature K = <• _ i\i 

on M^ having X[x,y) = (1,0), Y{x,y) = (0, y — x^) as orthonornial frame. The curvature has a 
crest passing through the tangency point (0, 0) 



open sets M~^ , M such that Z = dM^ = dM , f : -E|a/+ — )• TM~^ is an orientation-preserving 
isomorphism and f : E\j^,j- — )• TM~ is an orientation reversing-isomorphism. 

Remark 2 Throughout this section we deal with totally oriented ARS on compact surfaces (satis- 
fying assumption (HO), see RemarkU^. 



5.1 Gauss-Bonnet formulas 

The first global result for almost-Riemannian surfaces is a Gauss-Bonnet-type formula proved in 
[3]. Such theorem gives a first sight of how the relation between curvature and topology changes 
in the almost-Riemannian context. Indeed, not only the topology of the surface shows up in the 
Gauss-Bonnet formula, but also the way the singular set embeds in the surface. Moreover, a central 
role is played by the topology of the vector bundle, see Corollary [H 

Other generalizations of the Gauss-Bonnet formula can be found in [5] for contact sub-Riemannian 
manifolds and in [27^ [28] for pseudo-Riemannian manifolds. 



5.1.1 ARSs without tangency points 

The Gauss-Bonnet formula for a compact oriented Riemannian surface M states that 



[ KdA = 2TrxiM), 
Jm 



(10) 



where K is the Gaussian curvature and dA is the volume form associated with the metric. 

The first step to generalize (|10p to ARSs is to introduce an analogous of the volume form dA. 
To this aim, the idea is to take a non-degenerate volume form lo on E for the Euclidean structure 
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and push it forward on M using the morphism f. Since f is not an isomorphism on Z^ this does not 
give rise to a two form on the whole surface M but only on the set M\Z. Let dA^ = f^w denote 
such a form. Then dAg defines on M~^ the same orientation as the one chosen on M, the opposite 
one on M~ . 

By construction, dAs G ^'^{M\Z) and it diverges when approaching the singular set. Also, the 
Gaussian curvature diverges when approaching to Z. To overcome these issues and provide a way 
of integrating KdAg over M the idea in [3j is to integrate KdAg over the set of points e-far from 
the singular set and then let e go to zero. More precisely, let 

M, = {qeM\d{q,Z)>e}, (11) 

(where d is the almost-Riemannian distance). Then the Gaussian curvature is said to be integrable 
with respect to the ARS if the limit 

lim / KdAs (12) 

exists. In this case such limit is denoted by J^j KdAg. 

It turns out that if there are not tangency points then the limit in (J12p exists and can be 
calculated in terms of the topology of the ARS. 

Theorem 3 ([Sj) If there are no tangency points, then 

[ KdAs = 2^(x(M+) - x(M-)). (13) 

To explain the assumption of absence of tangency points, let us go through the main ideas in 
the proof of (jlSp . By definition, 

KdAs = I KdAs + I KdAs. 

A/e JKUr\M+ JMeHM- 

On the sets M^ n M^ we can apply the classical Gauss-Bonnet formula to get 

/ KdAs = 2TTxiMe n M+) - / kgda, 

where kg is the geodesic curvature and d{M^ n M"*") carries the orientation induced by M^ n M~^ . 
Taking account of orientations, similarly we obtain 

/ KdAs = -2Txx{Me n M~) + / kgda. 

Hence letting e go to zero, to show formula (113p it is sufficient to prove that the contribution at 
boundaries offset each other, i.e., 

lim ( / kgda- kgda] = 0. (14) 

This is the key point where the absence of tangency point is used. Indeed, under this assumption, 
the set dM^ is shown to be smooth and, moreover, the contributions at the boundaries have the 
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same order in e. In next section we see that this symmetry between the two sides of the singular 
set is lost in a neighborhood of a tangency point. 

Theorem [3] was generalized in [21j to surfaces with boundary. In this case, the authors define 
admissible domains as open bounded connected domains U <Z M whose boundary is the finite 
union of C^-smooth admissible curves and study the convergence of /r/^jv/ KdAg as e goes to zero. 
The generalized Gauss-Bonnet formula, proved through the above techniques, takes account of the 
boundary contributions and some other terms due to the intersections of dU with the singular set. 
When the intersection between dU and Z is C^-smooth, the formula simplifies to the following one. 

Theorem 4 ([21]) If there are no tangency points and U is an admissible domain such that dU 
is piecewise C" -smooth and dU Z is C" -smooth, then 



'U JdU 

In the previous statement we have 



/ KdAs + [ kgda = 27r(x(C/+) - xiU')). 

Ju JdU 



KdAs = lim / KdAg, 



/ KdAs = lim ( / kgdAg - / kgdA, 

JdU ^~^° \Jdundu+ Jaundur 

where U has the orientation induced as a domain of M and dU is oriented as boundary of U; 
U^ denotes the set U n M^ n M^ and kg denotes the Gaussian curvature and da the arclength 
parameteiiJ. 

5.1.2 ARSs with tangency points 

Let us present a further generalization of Theorem [3] allowing the presence of tangency point. 

In this case, two main issues are to be considered when studying the convergence (|12p . First, 
around a tangency point the boundary of the domain of integration is not smooth but has two 
corners (one on each side of Z, see Figure [6]). Second, the Gaussian curvature diverges in a more 
complicated way than at Grushin points. Indeed, while at Grushin points K diverges to — oo, at 
a tangency point K diverges at +oo on some directions, and it diverges at — oo at some other 
ones (see Figures El [9]) . Together with the interaction between different orders in the asymptotic 
expansion of the almost-Riemannian distance, these remarks possibly explain why the limit (J12p 
is not shown to converge. Indeed, the compensation between boundary terms ()14p appears not to 
happen around tangency points. This has been supported by numerical simulations in [HI Section 
5.1] for the ARS on M? having (1,0), (y — x^) as orthonormal frame. For this example, the limit 
(fHj) appears to diverge near (0,0) as c/e, c is a positive constant. 

To overcome the problem, a new notion of integrability of the curvature has been provided in 
[8] . The idea is to integrate the curvature not on the whole Mg , but on a subset of Me depending 
on two other parameters 61,62- This set is built by taking away from M^ a "rectangular" box for 
each tangency point, where 61, 62 are the dimensions of the box see Figure [TOl 



^Each C^-smooth piece of dU is an admissible curve parameterized by arclength 
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K <0 



K = 



Figure 9: Behavior of K for the ARS on M? whose orthonornial frame is (1,0), (0,y — x^): the 
plane R^ with the singular set (solid line) and the set of points at which K = (dashed line). See 
also Figure [8] for the graph of K 



7^2 (-^i) 




-y-S^i-Si) 



7^2 {5i) 



7-52 (5i) 



Figure 10: Let g be a tangency point. Take any smooth parameterized curve c :] — 1, 1[— )• M such 
that c(0) = q and c{s) ^ A(c(s)). For each s G (—1, 1), denote by 7^ the geodesic (parameterized 
by arclength) such that 7s(0) = c{s) and d{jsit),w{] — 1, 1[) = \t\ for each t sufficiently small. For 
61, 62 sufficiently small, the rectangle Bj ^ is the subset of M containing the tangency point q and 
having as boundary -fs2{[-^i,^i]) U 7[-52,52]('^i) U 7-<52([-'Ji, (^i]) U 7[_52,<52](-'5i) 



Let M^^Si,S2 = ^e \ UgeT^S s ■ Then K is said to be 3-scale integrable with respect to the 
ARS if the limit 



lim lim lim 

Si^O 62^0 e^O J j^j 



KdA, 



(15) 



exists. In this case such limit is denoted by J^i^j KdAg. 

Note that when there are no tangency points, the limit defined in (J15p clearly coincides with 
the one in (J12p . Besides, the order in which the limits are taken in (|15p is important. Indeed, if 
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the order is permuted, then the result given in Theorem [5] does not hold. 

To obtain a Gauss-Bonnet formula when tangency points are present, we need the notion of 
contribution at tangency points. 

Definition 4 Let q he a tangency point. Orient Z as the boundary of M~^ (see Figure [Tl]) . We 
define the contribution at q as Tq = 1, respectively Tq = —1, if the distribution is rotating counter- 
clockwise, respectively clockwise, along Z at q, see Figure [711 




M+ ^ \ \ / I M+ 




Figure 11: Tangency points with opposite contributions 

The following result generalizes Theorem [3] to ARSs with tangency points. 
Theorem 5 ([8j) The Gaussian curvature is 3-scale integrable and 



^ KdAs = 2^ x(M+) - X{M-) + 5] rj , 



(16) 



where T is the set of tangency points of S. 



Notice that the construction of M^^^^^^j depends on the choice of a manifold transversal to Z 
at each tangency point and on its parameterization. A canonical choice of this manifold has been 
provided in |18l Theorem 2], implying that M^^g-^^g^ can be constructed in a intrinsic way. This is 
also suggested by the fact that [|i^^ KdAg equals a quantity that is intrinsically associated with the 
structure, see formula (fTHl) . 

Proof of Theorem [5] (sketch). Fix 5i and 82 in such a way that the rectangles Bg g are pairwise 
disjoint and Z n dB^ ^ C [—5i, 5i] x {82}, for every q G T. By construction, dB^ ^ is admissible 
and has finite length for every q €z T- Hence we can take M^ \ U^eT ^5 5 as C/ in Theorem HI As 
a consequence, we have 



/ KdAs + J2[ kgdas = 27r{x{M+\\jBlg^)-x{M-\\jBlgJ) 

= 27r(x(M+)-x(M-)). (17) 



If we prove that, for a fixed q £ T, 



lim lim / k„das = — 27rr„, (18) 
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then we directly obtain ([T6]). To deduce ([T8|) consider the local representation {F3) and assume 
that {{x^y) I y — x^ip{x) < 0} C Af+, the proof for the case {{x,y) \ y — x'^'4>{x) < 0} C M~ being 
analogous. On one hand, one can check that Tq = 1. On the other hand, the geodesic curvature 
along [— (5i,(5i] x {62} and along [— 5i,5i] x {—62} is zero, the two segments being the support of 
geodesies. Hence 

/ kgdas = I kgdas + / kgdag + T^ Oj 

JdBl^,^ J{5j}x[-52,S2] J{-5^}x[-52,S2] p^ 

where the last term is the sum of the values of the angles of the box and is equal to —2tt. Indeed, 
because of the diagonal form of the metric with respect to the chosen coordinates, each angle has 
value — ^. The first two terms are well defined and tend to zero when 62 tends to zero. Hence 

lim lim / kgdus = — 27r = —2-kt„. 

Si^O 62-^0 J QB'' ^ 



5.2 A topological classification of ARSs 

In this section we present a result that provides a relation among the topology of the almost- 
Riemannian surface and the Euler number of the vector bundle associated with the structure. 

Let us recall the notion of Euler number. Given an oriented vector bundle of rank 2 over 
a compact connected oriented surface M, the Euler number of E, denoted by e{E), is the self- 
intersection number of M in E, where M is identified with the zero section. To compute e{E), 
consider a smooth section a : M ^ E transverse to the zero section. Then, by definition, 

e{E)= Yl ^(P'^)' 

p\a{p)=0 

where i{p,a) = 1, respectively —1, if dpa : TpM — )• T^^^p-^E preserves, respectively reverses, the 
orientation. Notice that if we reverse the orientation on M or on S then e{E) changes sign. 
Hence, the Euler number of an orientable vector bundle E is defined up to a sign, depending 
on the orientations of both E and M. Since reversing the orientation on M also reverses the 
orientation of TM, the Euler number of TM is defined unambiguously and is equal to xi^)^ the 
Euler characteristic of M. We refer the reader to [24] for a more detailed discussion of this subject. 

Theorem 6 ([8j) Under the assumptions of Remark\^ there holds 

X{M+) - x(M-) + 5^ T, = e{E). (19) 

The strategy to prove (|19|) is based on the construction of a section a oi E having only isolated 
zeros {pi, . . . ,Pm\ and such that 

m 

Y, i{Pi,'y) = x{M+) - x(M-) + Y, ^1- 
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To this aim, the key point is to construct cr in a tubular neighborhood of Z in such a way that it 
vanishes only at tangency points (see [8, Lemma 1]) and at each g S T the relation i{q,cr) = Tq 
holds (see [8i Lemma 2]). Once this is done, it is sufficient to extend o" in a smooth way to the 
whole surface. By a transversality argument, this can be done by introducing only a finite number 
of isolated zeros whose index sum can be calculated using Hopf 's Index Formula. 

Theorem [6] has several implications. First it classifies the ARS with respect to the associated 
vector bundle. Indeed, the Euler number represents the only topological invariant of an oriented 
rank-2 vector bundle over a compact oriented surface, i.e., it identifies the vector bundle. In 
particular, as a direct consequence we get that an ARS is trivializable, i.e., E is isomorphic to the 
trivial bundle, if and only if 

x(m+)-x(m-) + J;t, = o. 

This generalizes and provides the converse result of [HI Lemma 5] stating that if tangency point are 
absent, i.e., 7" = 0, and the structure is trivializable then x(M+) — x{^~) = 0. An alternative 
proof of the fact that the latter condition is sufficient for the structure to be trivializable can be 
found in [7j. 

Moreover, Theorem [6] allows to rewrite the Gauss-Bonnet formula of Theorem [5] as follows. 

Corollary 1 For any totally oriented ARS S on a compact surface M , the Gaussian curvature is 
3-scale integrable and 

m KdAs = 27re{E). (20) 

Tm 

Corollary [1] encloses previous Gauss-Bonnet formulas: for Riemannian structures, where f is an 
isomorphism and |e(S)| = e{TM) = x{^) (formula (llOp ) and for ARSs without tangency points, 
where T = (formula (fT3]) ). 

When considering Riemannian structures formula (j20p says that the total curvature is zero if 
and only xi^) = 0) that is, M is diffeomorphic to the torus. Instead, in the almost-Riemannian 
context formula (I20p implies that if the total curvature is zero then the vector bundle E is trivial. 
As a consequence, there exist ARSs on surfaces of positive genus having zero total curvature, see 
for example Section! 



5.3 Lipschitz equivalence of ARSs 

This section is devoted to the description of how the presence of the singular set and, in particular, 
of tangency points affect the distance associated with the ARS. Namely, we focus our attention on 
the problem of Lipschitz equivalence among different almost-Riemannian distances. 

A Lipschitz equivalence is a diffeomorphism if : Mi — t- M2 which is bi-Lipschitz as a map from 
the metric space {Mi,di) to {M2,d2), where di is an almost-Riemannian distance on the surface 
Mi associated with an ARS {Ei,fi, (•, •)j) on Mj. Recall that bi-Lipschitz means that there exists a 
constant C > 1 such that 

— d2{(p{q),ip{p)) < di{q,p) < C d2{ip{q) , ip{p)) , ^ q,p ^ Mi. 

In the Riemannian case, all distances on diffeomorphic compact oriented surfaces are Lipschitz 
equivalent. In other words the Lipschitz classification of Riemannian distances on compact oriented 
surfaces coincides with the differential one. 

22 



In the almost-Riemannian case the Lipschitz classification is finer. Clearly this is due to the 
presence of a singular set and mainly it is due to how the singular set splits the surface together 
with the location of tangency points with their contributions. 



5.3.1 Graph of a totally oriented ARS 

It turns out that all the information needed to identify the Lipschitz equivalence class of an almost- 
Riemannian distance can be encoded in a labelled graph that is naturally associated with the 
structure. 

The vertices of such graph correspond to connected components oi M \Z and the edges corre- 
spond to connected components of Z. The edge corresponding to a connected component W oi Z 
joins the two vertices corresponding to the connected components oi M \Z adjacent to W. Every 
vertex v corresponding to a component My is labelled with a pair of integers (sign(f ), x(^^)), where 
sign(t') takes into account of the orientation of M„ (sign(t;) = ±1 if My C M^) and x(^) is the 
Euler characteristic of M^. Every edge e corresponding to a component W C ^ is labelled with the 
ordered sequence of signs (modulo cyclic permutations) given by the contributions at the tangency 
points belonging to W where the order is fixed by walking along W oriented as the boundary of 
M+. In Figure [12] we illustrate the algorythm to build the labelled graph associated with an ARS. 



{E3- 






Figure 12: Algorythm to build the graph 

An example of ARS and associated graph is portraiteqj in Figure [T3l In this case we consider 
an ARS on the compact oriented surface of genus 2 where Z is the union of two circles placed as 
in Figure [T3b and where there is only a tangency point of positive contribution on one of them. 
Then the graph (Figure [T3b) must have 3 vertices and 2 edges. The label on vertices are easily 
computed. As concerns labels on edges there is only a tangency point (with positive contribution) 
on a connected component of Z, whence only one edge carries a label (1), while the other edge is 
labelled with (0). 



*Note that in Figure [TH] we represent only the surface, the singular set and the tangency points. Here we are not 
interested in the explicit expression of the morphism f nor on the vector bundle E. 
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Remark 3 Adding the quantity ^^ sign(t;)x('u) to the sum of all the entries in the label of edges, 
we obtain x{M~^) — x(M~) + X^„g7-Tg, which equals e{E), by Theorem\^ Also, YlvX{v) = x{^)- 
Hence, once the labelled graph associated with an ARS is given, one recovers directly the vector 
bundle and the surface. 

Moreover, the labelled graph associated with an ARS depends on the orientation fixed on E. 
More precisely, choosing on E the opposite orientation produces the following changes in the labels 
of the graph. On each vertex the first entry of the label changes sign. On each edge not only each 
entry of the tuple changes sign but also the tuple is reversed in order]^ 



r = 





-1,-2 



Figure 13: Example of ARS and associated graph 

We say that two labelled graphs associated with two totally oriented ARSs are equivalent if 
either they are equal or after possibly changing the orientation on one vector bundle they are equal. 
This notion of equivalence is motivated by this straightforward remark: changing the orientation 
of the vector bundle does not affect the almost-Riemannian distance. Moreover, when two labelled 
graphs are equivalent, by Remark [3] the associated vector bundles are isomorphic and the underlying 
surfaces are diffeomorphic. 

5.3.2 A classification result 

The following result classifies totally oriented ARSs with respect to Lipschitz equivalence. 

Theorem 7 ( |19] ) Two totally oriented almost-Riemannian structures defined on compact sur- 
faces are Lipschitz equivalent if and only if they have equivalent graphs. 

Theorem [7] implies that the Lipschitz equivalence class of an almost-Riemannian distance on 
a surface depends only on how the singular set is embedded in the surface and how the tangency 
points with their contribution are located. 

In the Riemannian contexl|^, the Lipschitz equivalence class of distances on a given surface 
is unique or, equivalently, it does not depend on the bilinear form on the tangent bundle. In a 
similar way, by Theorem [7] we deduce that the Lipschitz equivalence between two distances does 



'For example if the label of an edge is ( — 1, 1, — 1, — 1), changing the orientation on E the new label becomes 
(1,1,-1,1). 

^"if the structure is Riemannian, the associated labelled graph consists of a unique vertex and has no edges. The 
label on the vertex is {5,x{M)), where 5 = 1 if f preserves the orientation, —1 otherwise. 
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not depend on the bilinear form G defined on A (see Section [2|) but only on the submodule A 
itself. This is highlightened by the fact that the graph itself depends only on A. In terms of the 
triple (-E, f, (•, •)) this translates to the fact that the labelled graph does not depend on the chosen 
Euclidean structure (•,•), but only on the morphism f. As a consequence, in general Lipschitz 
equivalence does not imply isometry. 

The main tool in the proof of Theorem[7]is a local classification of ARS by Lispchitz equivalence. 
Indeed, using the Ball-Box Theorem (see \1Q\ Corollary 7.35]) it is not hard to show that in a 
neighborhood of a given point two ARSs are Lipschitz equivalent if and only if the point is of the 
same type (Riemannian, Grushin, tangency) for both structures. Then, to glue the information 
on different neighborhoods one needs the topological classification of surfaces, that is, one uses 
the information carried by labels on vertices. To appreciate the role of contributions of tangency 
points, recall that in a neighborhood of a tangency point the asymptotic of the almost-Riemannian 



distance from Z is different from the two sides of Z, see Section 14.21 The contribution Tq carries 
the information about which side is the one with the fastest rate of convergence of Mg (see (llip ) 
to Z. 



6 Conclusions 

In Table 1 we compair Riemannian and almost-Riemannian geometry of surfaces summing up the 
main aspects presented in this paper. 



RIEMANNIAN 


ALMOST-RIEMANNIAN 


X, Y linearly independent 


X, Y Lie bracket generating 


at each point spheres of small radii are smooth 


at points in Z spheres are never smooth 


K < then no conjugate points 


both Riemannian points and points in Z have conjugate locus 


the cut locus from a point is separeted from the point 


the cut locus from a point in Z accumulates at the point 


Lipschitz equivalence = differential equivalence 


Lipschitz equivalence C differential equivalence 


Jj^KdA = 27rx{M) 


#^ KdAs = 2ne{E) 


trivializable (3 a global orthonormal frame) =^ torus 


trivializable =► torus or Z y^ 9 



Table 1 : Comparison between Riemannian and almost-Riemannian geometry 



As one can infer from the present analysis, the most interesting points of ARSs are tangency 
points. Even though some contributions have been done |131 I18j . tangency points are far to be 
deeply understood. 

An open question arisen in the proof of Theorem [3] is the convergence or the divergence of the 
integral of the geodesic curvature on the boundary of a tubular neighborhood of the singular set, 
close to a tangency point. To address this issue one needs to understand whether the domain of 
integration or the form to be integrated need to be reconsidered. 

The analysis of the Laplace-Beltrami operator in presence of tangency points has not been 
considered yet although in [20] the authors conjecture some properties, based on the case with only 
Grushin points. 
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Another interesting problem is whether it is possible to associate with an ARS a canonical linear 
connection on E compatible with the Euclidean structure. To this aim one needs to use properties 
of the morphism f , as in general the answer is negative. This could be a step forward towards an 
intrinsic notion of integration of the Gaussian curvature on the surface, as one could study the 
curvature on the vector bundle. Also, focusing on the vector bundle rather than on the (tangent 
bundle to the) surface is a starting point towards higher dimensions generalizations. 
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